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ABSTRACT 


The thermal performance characteristics of an electronics cooling Liquid 
Flow Through Module (FTM) were experimentally investigated. Different sets of 
experiments were conducted for each side of the FTM. A synthetic dielectric 
polyalphaolefm type coolant, Brayco Micronic 889, was used. Six etched foil type 
heaters were attached to one side of the FTM over the fluid flow path while three 
heaters were attached to the other side of the module. Inlet and outlet fluid 
temperatures as well as surface temperature data were acquired from both sides of 
the module for several different flow rate and power setting combinations to 
quantify the effectiveness of the FTM. Correlations, in terms of Reynolds and 
Stanton munbers were formulated according to the data for both sides of the 


module. 
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INTRODUCTION 


A. OVERVIEW 

Since the development of the first integrated circuits (IC), there has been a steady 
increase in heat dissipation levels of electronic devices. This increase has made the role 
of heat transfer and thermal design more important than ever. The necessity for advanced 
techniques in electronic cooling started with large-scale integration (LSI) technologies of 
the 1970's in which 1000 junctions per chip were achieved. This continued with very 
large-scale integration (VLSI) technologies of the 1980's in which 100,000 junctions per 
chip were achieved [5]. The magnitude of this thermal stress problem becomes more 
apparent when microelectronics trends are examined. Microprocessor clock frequencies 
have been doubling every five years and as clock speeds increased, almost linearly, so did 
heat dissipation levels. Today’s technology has achieved 75-100 MHz system clock speed 
or even higher. Therefore 400-500 MHz clock speed could be expected towards the end 
of this century. Thus a 500 MHz electronic device will dissipate 5 times the heat 
dissipated at 100 MHz [10]. 

These trends indicate that thermal stress management will always become a major 
problem in electronic packaging. Higher thermal stresses will result higher failure rates, 
shorter operational lifetimes and higher logistics and maintenance costs. 

There are several advanced thermal control techniques to improve cooling 
efficiency. Since the improvements in reliability and packaging density of electronic 
devices are directly related to advances in thermal design, the choice of thermal control 
technology has large effects on reliability and cost. Insufficient attention to thermal 
control may easily result in device failure. On the other hand thermal overdesign may 
result very high product costs. By taking care of these two criteria: reliability and cost. 



many new thermal control techniques involving liquid cooling, have been proposed. 
These new liquid cooling techniques can be grouped as follows; 

a. Jet Impingement and Nucleate Boiling Cooling 

b. MicroChannel Cooling 

These techniques will be discussed briefly in the following sections. 

B. AIMS OF THERMAL CONTROL 

The primary aim of thermal control is to prevent catastrophic failure which can 
be defined as immediate, thermally induced, total loss of electronic function. Catastrophic 
failure may result from melting or even vaporization of a component in microelectronic 
devices and in particular may also result from plastic flow of printed circuits. 

Catastrophic failure is generally dependent upon the local temperature field as 
well as the operating history. Therefore it is difficult to precisely determine where such 
failure can be expected to occur. However, maximum allowable temperatures can be 
found in manufacturer’s catalogs by which the possible range of component variations can 
be taken into account. Likewise, the possible range of environmental variations is another 
prerequisite for developing designs of thermal control systems [5]. 

C. ADVANCED LIQUID COOLING TECHNIQUES 

Natural Convection and radiation cooling, relying upon the ambient air, is 
insufficient to deal with the rapidly increasing heat fluxes. Figure -1 shows the 
relationship of the tremendous increase in Microelectronics heat fluxes to other 
technologies [5]. 

On the other hand, conventional forced convection air cooling can no longer 
satisfy heat removal requirements. As shown in Figure -2 , at lower air velocities, typical 
heat flux limits, for packages cooled by conventional forced convection, are less than 1 
W/cmMl2]. 
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Figure 1: Perspective on Microelectronic 
Heat Fluxes 



Figure 1. Perspective of Increase in Microelectronic Heat Fluxes [5], 

In the near future, technology will be able to produce very-large scale integrated 
(VLSI) chips dissipating up to 200 W/cm^ [7], To achieve a board level heat flux that can 
be cooled with conventional forced convection, package surface area must exceed total 
chip surface area by more than 25 times. For example, a forced convection cooled 
package containing 10 chips having 1 cm^ area each, would require 250 cm^ of surface 
area. Therefore it is obvious that conventional forced convection will not be applicable 
for today's VLSI electronic systems. However conventional air cooling is an important 
milestone which leads to the development of alternative methods of cooling. Temperature 
differences attainable as a function of heat flux for various heat transfer modes, and 
various coolant fluids is shown in Figure -3 [5]. 
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Figure 2. Package Heat Flux in Air Cooled Packages [12], 

1. Jet Impingement and Nucleate Boiling Cooling 

Alternative approaches to achieving high heat flux removal are to apply liquid - 
jet impingement, with or without boiling, or direct immersion cooling, so called "pool 
boiling", to the backside of silicon chips. Convective heat transfer associated with gas jets 
has been a research topic since the mid- 1950s, but jet impingement cooling with 
dielectric liquids is a recent improvement. In this cooling method there are two modes. 
Liquid may flow to the heater through a vapor / gas space or through a liquid - filled 
space. The first is called a "free-jet", and the latter is called a "submerged jet" [11]. High 
convective heat transfer coefficients can be achieved in this latter cooling mode. Thus 
"submerged jet" impingement cooling is one of the best alternatives in direct liquid 
cooling of electronic components. 
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Figure 3. Temperature Differences Attainable as a Function of Heat Flux 

for Various Heat Transfer Modes and Various Coolant Fluids [5], 

2. MicroChannel Cooling 

One of the best known cooling methods is the microchannel approach. In this 
technique, heat fluxes approaching 800 W/ cm^ with a temperature rise of TO^C have been 
achieved by using water as the working fluid [12]. These microchannels mentioned here, 
can be formed by plain fins, wavy fins, pin fins or ojfset strip fins. Offset strip fin type is 
the one which is widely used in microchannel designs. 

In this study, a Standard Electronic Module (SEM) which is manufactured by 
Lockhart Industries, has been investigated. Also in this Flow Through Module (FTM), 
which is a format- E size, offset strip fins have been used to increase the heat transfer 
coefficient by interrupting the flow so that the fluid boundary layer can never become too 
thick to reduce convective heat transfer. The flow is also more uniform than straight fins 
since it is a developing laminar flow through the microchannels formed by offset strip 
fins [8]. Typical rectangular offset plate fin exchanger and general notations of offset 
plate fins are shown in Figure - 4 [7]. 
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Figure 4. Typical Rectangular Offset Plate Fin Heat Exchanger and General 
Notation Used for the Dimensions of an Offset Plate Fin [7], 

Many experimental investigations have been conducted to improve the convective 
heat transfer coefficients by using fins. These studies dealt mostly with gases. Therefore 
the usefulness of these data is very limited when applied to liquid coolants, because 
Prandtl number which is a fluid property, has a significant effect on heat transfer 
coefficient. These tests have been conducted by Lockhart, for different fin pitches and fin 
heights, for Reynolds numbers varying between 30 to 2700 and Prandtl numbers from 0.7 
6 






to 70. The test results indicated that the Colburn j and the flow friction factors, are 
different for selected types of fluids within the same Reynolds number range [8], 
Therefore Kays' and London's air test data, [9] is not applicable to fluids which are 
outside the gas Prandtl number range, since the test results obtained by Lockhart confirm 
the Prandtl number effects. 

On the other hand there may be some disadvantages associated with liquid 
cooling, such as additional assembly relating to coolant circulation and refrigeration 
systems. Also the possibility of leakage may cause damage to the electronic components. 
That is why dielectric coolants are used in liquid cooling systems for electronics. In this 
study Brayco Micronic 889, manufactured by Castrol, has been used. This liquid is a 
polyalphaolefin (PAO) made for use in closed loop systems. Properties of Brayco 
Micronic 889 are given in Appendix-C. 

3. Liquid Flow-Through Modules 

In this study the effectiveness of a liquid Flow-Through Module (FTM) 
manufactured by Lockhart Industries, using PAO, is investigated. Flow-Through modules 
are manufactured by a technique so called "vacuum brazing" which facilitates the 
construction of microchannels in a very thin module. The module tested, shown in 
Figure-5, has a thickness of 0.100 inch , with 0.020 inch Aluminum 6061 faceplates. 

The fins within the module are made from an alloy. Aluminum 3003 and have 
0.06 inch height, 0.125 inch length and 0.006 inch thickness. The fin density is twenty 
fins per inch. [2] A schematic drawing of the fins is shown in Figure -6. 

Both sides of the FTM, shown in Figure -5, can be used for cooling electronic 
components. Electronic components placed on the surface of the FTM are cooled by the 
coolant flowing through the module. The thin structure of the module is very helpful for 
stacking purposes. Several modules can be housed in a small enclosure. For example, in 
a study performed by AT&T Bell Laboratories about Liquid Flow-Through cooling for 
avionics applications, 52 FTM's have been housed in a small enclosure. Both the module 



and the system have demonstrated very high cooling efficiency and good reliability under 
shock and vibration testing [10]. The FTM tested in this study, shown in Figure 7, 
consists of three main flow paths of microchannels formed by offset strip fins. The flow 
pattern of the coolant is shown in Figure -7 [2]. The FTM is currently under evaluation 
by the military for naval and avionics applications [1]. 

D. OBJECTIVES 

This study experimentally investigated the capacity for heat removal of the 
SEM-E sized Flow-Through Module, using a synthetic dielectric polyalphaolefin coolant, 
so called Brayco Micronic 889 , manufactured by Castrol. Experimental investigation was 
conducted for both sides of the FTM by utilizing several different flow rate and power 
setting combinations. Temperature data were collected from both sides of the FTM in 
order to quantify the effectiveness of the module. Correlations, in terms of Reynolds and 
Stanton numbers were formulated. According to the temperature data collected from the 
module inlet and outlet, heat dissipation capabilities of the FTM were quantified. 





Figure 5. Top View of the Flow-Through Module. 

















Figure 6. Dimensions of the Offset Strip Fins in the Flow-Through Module [2] 






FLOW DIRECTION 



Module Thickness : 0.1 inch 
Faceplate Thickness : 0.02 inch 
Faceplate Matenal : 6061 Aluminum 

Figure 7. Fluid Flow Path within the FTM [2], 
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11. EXPERIMENTAL APPARATUS 


A. SUPPORTING SYSTEM ASSEMBLY 

The schematic diagram of the experimental setup, as shown in Figure -9 can be 
examined in two parts. The first part, which is called " Supporting System Assembly" 
consists of: Fluid Circulation, Power Distribution and Data Acquisition sub-systems. Each 
of these sub-systems and their components are as follows. 

1. Fluid Circulation System 

This section consisted of refrigeration bath, pump, variable speed control unit and 
flowmeter. The refrigerated circulating bath which provided the heat sink for the coolant 
was an Endocal RTE- 5. The pump and the drive speed control unit, manufactured by 
Masterflex could provide a maximum volumetric flow rate of 1 liter per minute. Tbe 
coolant was pumped through flexible Tygon tubing. A flowmeter manufactured by Fischer 
and Porter was used to determine the flow rate. The Tygon tubing was connected to the 
Flow-Through Module through brass tubing, attached to the inlet and outlet ports of the 
module. The brass tubes were inserted into Tygon tubing and clamps were affixed to 
tighten the connections. Upon exiting the module the coolant returned to the refrigerated 
bath. 


2. Power Distribution System 

A Hewlett Packard 0-20V, O-lOA DC power supply was used to supply power to 
the heaters. Six heaters were affixed to the side-A and three heaters were affixed to the 
side-B, by using thermally conductive epoxy adhesive. Omega Bond 101. The heaters, 
manufactured by Minco, all wired in parallel, were 7.6 mm (0.3 in) by 38.1 mm (1.5 in) , 
kapton etched foil heaters. 
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itai Setup. 







The six heaters on side-A, as shown in Figure -10, were placed 5 cm (2 in) apart 
along the module flow paths by affixing two heaters on each flow path. The three heaters 
on side-B, as shown in Figure -11, were placed on the centerline of the module, the rows 
of heaters were 1.25 cm (0.5 in) apart from each other. Each heater had an effective 
heated area of 1.94 cm^ (0.3 in^ and a resistance of 20.7 ohms ± 10 %. 

In order to measure the current passing through each heater, thermally stable 
precision resistors were connected in series with each heater. Prior to connecting, each 
resistor was calibrated by using a Rosemount 920 A Commutating Bridge to determine its 
resistance precisely. Six of the resistors, connected in series with the six heaters on side-A 
were 0.487 ohms ± 1 % and the other three resistors connected in series with the three 
heaters on side-B were 0.500 ohms ± 1 %. 

A Hewlett Packard 44701A digital voltmeter which was a part of the HP 3852A 
Data Acquisition Control Unit, was used to determine the voltage drops across the heaters 
and the precision resistors. All of the 18 measurements, voltage drops across the heaters 
and precision resistors were measured by utilizing an HP 44705A 20 channel multiplexer. 

3. Data Acquisition System 

A Hewlett Packard 3852A Data Acquisition / Control Unit was used to measure 
the voltages and temperatures. In order to be able communicate with the data Acquisition 
/ Control Unit through a computer, an interface board, GPIB-PCIIA IEEE-488, 
manufactured by National Instruments had been installed in a Personal Computer. A 
program written in QBASIC code, attached as Appendix-D, was used to acquire 
temperature and voltage data. An HP 44708A 20 channel relay type multiplexer was used 
for temperature measurements. The experiments were conducted first for side-A and then 
repeated for side-B, since all the 20 chatmels of the HP44708A were utilized during the 
experiments. Thermocouple locations for both sides of the FTM were shown, 
schematically in Figures 12-a and 12-b. 

Fifteen of the thermocouples were attached onto the module itself and two 
thermocouples were placed at the inlet and the outlet of the module. Two thermocouples 
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Figure 11. Heater Locations on the rTM for Side - B, 
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Fifteen of the theraiocouples were attached onto the module itself and two 
thermocouples were placed at the inlet and the outlet of the module. Two thermocouples 
were placed in the insulation layers to estimate the heat loss. An additional thermocouple 
was used either to measure ambient room temperature or coolant temperature in the 
refrigerated bath. Thermocouple extensions were fed into channel numbers 0 through 19 
of the HP44708A in the data acquisition system. 

Voltage drops across the heaters were measured by the charmels 0 through 8 of 
the HP44705A multiplexer, corresponded to heaters A through I respectively. On the 
other hand , precision resistor voltages were measured by utilizing the channels 9 through 
17 of the HP44705A mutiplexer. 

B. TEST SECTION ASSEMBLY 

1. Housing of the Flow-Through Module 

The Flow-Through Module was supported in a horizontal position by a plexiglass 
enclosure. The FTM was insulated on all sides, inside the plexiglass enclosure, by using 
0.75 inch (1.91 cm) thick, black foam type insulation. Two layers of insulation were 
placed on top and bottom sides of the FTM and one layer of insulation was placed along 
the edges of the module as well. The front side of the enclosure was open for wiring and 
tubing connections [1]. 

2. Implementation of Thermocouples 

Thermocouples affixed onto the surface of the module were, T-type 
Copper-Constantan, flat, 5 mil thick, foil type, cement-on thermocouples. They could be 
placed very close to the heaters. Copper leads were soldered to the thermocouple 
extension wire and the constantan leads were resistance welded to the extension wire. 
These connections were wrapped using kapton tape to keep them separate. 
Thermocouples were affixed by using the same kapton tape. 



T- type thermocouple probes were used to measure the coolant temperature at the 
inlet and outlet of the module. These thermocouples were 8 mil thick and are welded to 
form a "bead" junction. 

Thermocouples used to measure the temperatures at the top and bottom sides of 
the insulation were also 8 mil thick are welded T- type copper constantan thermocouples. 
They were affixed onto the insulation by using electrical tape. 

Thermocouples used at the inlet and outlet of the module for coolant temperature 
measurements, were calibrated by using a Rosemount 920A commutating bridge and a 
constant temperature ethylene glycol bath. 


21 



22 



III. EXPERIMENTAL PROCEDURES 


A. CALCULATION OF COOLANT THERMOPHYSICAL PROPERTIES 

Temperatures at the inlet an outlet of the module were measured 30 times during 
each run in a 20-second time interval. The average values of both inlet and outlet 
temperatures were computed and the mean of these two values was utilized to calculate 
the thermophysical properties of the coolant, Brayco Micronic 889. 


1. Density (g/ml) 

Densities of the coolant at various temperatures provided by the manufacturer 
(Castrol) were curve fitted and the polynomial obtained by this process was used to 
calculate the density at any temperature. Data Supplied by Castrol is tabulated in Table 1. 


TEMPERATURE ( ° C ) 

DENSITY (g / ml) 

0 ° C (32 ^ F) 

0.811 

20°C(68°F) 

0.794 

40 “C (104° F) 

0.777 

100‘’C(212°F) 

0.723 

160 “C (320 ° F) 

0.661 


Table 1. Density Data of the Coolant, Brayco Micronic 889 [3] . 

The curve obtained is shown in Figure 13. The second order polynomial the 
curve fit is as follows: 

p = -8.8216* 10-''* 7lvg-7.9251 * +0.81055 g/ml (1) 

where T,,^ is the average temperature differential of the coolant between inlet and outlet 
of the FTM in degrees Celsius. 
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Figure 13. Density Curve of the Coolant Brayco Micronic 889. 

2. SpeciHc Heat (J / g "C) 

In the previous study [1], certain inconsistencies in the energy balancearose. After 
thorough investigation it was hypothesized that the values for specific heat of the PAO 
were not accurate. To obtain accurate values of the specific heat, samples of the PAO, 
were sent to Delsen Testing Laboratories lnc.,for testing and measurement of the specific 
heat. A second order polynomial curve fit to this data is given as follows: 

Cp = 3.5284 ♦ 10-^ ♦ 71,^ + 1.0466 * lO'^ * r„vg + 0.46955 cal/gT (3) 

A comparison of the Delsen Laboratories data with the PAO manufacturer's data is 
given in Table 2 and Figure 14. 
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TEMPERATLIRE 

SPECIFIC HEAT 
(CASTROL, cal/g ° C) 

TEMPERATURE 

SPECIFIC HEAT 
(LAB'S DATA, cal/g ° C) 

-18 °C (O-’F) 

0.49 

I0°C(50°F) 

0.48 

lO-CCSO-'F) 

0.52 

38°C(I00°F) 

0.51 

38° C(100°F) 

0.54 

66°C(I51“F) 

0.54 

93°C(200°F) 

0.58 , 93 ° C (200 ° F) 

0.57 


Table 2 . Specific Heat Values of the Coolant (by Castrol and Delsen Laboratories). 
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Figure 14. Specific Heat Comparison of Brayco Micronic 889. 
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3. Kinematic Viscosity (mVs) 

Manufacturer's data was utilized to acquire a relationship for kinematic viscosity of 
the coolant in the temperature range of the experiments. Data supplied by Castrol is 
tabulated in Table 3 and plotted in Figure 15. With some uncertainty kinematic viscosity 
was assumed to be changing linearly in this temperature range. According to the data 
shown in Appendix -C, the relationship for the kinematic viscosity of Brayco Micronic 889 
was determined as follows; 



Table 3. Kinematic Viscosity Data of Brayco Micronic 889 [3]. 



Figure 15. Kinematic Viscosity Curve of Brayco Micronic 889. 
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B. FLOWMETER CALIBRATION PROCEDURE 

Flowmeter calibration was performed by measuring the time intervals for 1 and 
two liters of coolant to fill a scaled container for 10 percent increments of the flow rates 
from 10% to 100%. The data obtained are tabulated in Table 4. The average of two of 
the resulting measurements was taken and the data obtained was curve fitted. The 
calibration curve of the flowmeter for the coolant, Brayco Micronic 889 is shown in 
Figure 16. The fourth order polynomial curve fit is given below: 

S = 2.7178 xlO-’a''-6.02 xl0-V+0.422a^+4.1837a-8.07 ml/min (5) 
where a is the flowmeter reading in percentage. 



Figure 16. Flowmeter Calibration Curve for the Coolant, Brayco Micronic 889. 
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C. COLLECTION AND ANALYSIS OF THE DATA 

In this study two different sets of experiments were conducted: one for Side- A 
and one for Side- B of the Flow-Through Module. The naming of the sides of the module 
as A and B was arbitrary For the first set of experiments, thermocouples were affixed to 
Side- A where 6 heaters were located. Nine different power settings from 20 W to 100 W 
in 10 Watt increments and seven different flow rates from 350 ml/min to 950 ml/min in 
100 ml/min increments were tested. 

During the second set of experiments, which were for Side- B where 3 heaters 
were located, 5 different power settings from 20 W to 100 W in 20 Watt increments 
were tested. The flow rates were the same as those used for the first set of experiments. 
Temperature data were observed on the monitor of the PC to determine when steady state 
had been achieved. Steady state was assumed to be achieved when the average surface 
temperature and the temperature difference of the coolant between the thermocouples 
located at the inlet and the outlet of the module varied by no more than 0.02 “C over a 
5-rainute period. Voltage measurements, thermocouple readings, and a series of 30 
coolant inlet and outlet temperature readings were taken over a period of 20 seconds and 
the data was recorded. Depending upon the power setting and flow rate combination, 
steady state was achieved in 10 to 40 minutes. After saving the data to a file, the flow 
rate was adjusted to the next value to continue the procedure. 

It was generally observed that steady state was more rapidly achieved as the flow 
rate got higher. Flow rates were increased in 100 ml/min increments from 350 ml/min to 
950 ml/min. The flow began to pulse at lower flow rates. Therefore 350 ml/min was 
taken as the lower limit for the flowmeter settings. 

Power settings were adjusted using the HP power supply. The total amount of 
power supplied was monitored and adjusted until the desired value was achieved. Power 
settings were varied in the range from 20 W to 100 W by either 10 W or 20 Watt 
increments The diagram of the circuit formed by heaters and precision resistors was 
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shown in Figure-17, The variation of resistances could be 10 % in the heaters and 1 % in 
precision resistors according to the manufacturers' data. 



Figure 17. Schematic Drawing of the Electrical Circuit Formed by Heaters and Precision 
Resistors. 


1. Calculation of Total Power 

Data was analyzed by calculating the total power supplied by the heaters and 
comparing it with the sum of heat dissipated by the coolant and the heat lost through the 
insulation. 

Pt=qnu,d + qioss (6) 
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Current passing through each heater was calculated by utilizing the precision 
resistors, since each heater was eonnected in series as shown in Figure-17 

Vp 

= (7) 

Thus 9 different current values were obtained as follows; 

/*(0 = 4(0 = ^ i-0to8 (8) 

Then these current values used in the following equation to get the power 
supplied to each heater, 

F/,(/)= Ka(0x/a( 0 i = 0to8 (9) 

As a result the total power supplied to the module was calculated as follows; 

Watts (10) 

The heat flux generated by each heater in Watts/cm^ was calculated for each 
power setting by the following equation: 

= (Watts/cm= ) (11) 

Aeff y 

where A^j- = 1.94 cm^ Effective Heater area. 

The values of total power and heat flux per heater for each power setting were 
tabulated in Table 5. 
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Power Setting 
Side - A 

Power Setting 
Side - B 

P,«a,(W) 

Heat Flux per 
Heater (W/cm^) 

1 

1 

20 

1.15 

2 


30 

1.72 

3 

2 

40 

2.3 

4 


50 

2.86 

5 

3 

60 

3.44 

1 6 


70 

4.01 

7 

4 

80 

4.58 

8 

1 90 

5.15 

9 

5 1 100 

5.73 


Table 5 . Power Settings for Both Set of Experiments. 


2. Calculation of Heat Dissipation 

a. Average Temperature Differential 

The temperature of the coolant was measured at the inlet and outlet of the 
module for 30 times during each run of experiments. The average of both inlet and outlet 
temperatures were calculated by the program written in QBASIC for the data acquisition 
system as follows; 




( 12 ) 


b. Heat Dissipated by the Coolant 

The coolant dissipated almost all of the heat supplied by the heaters to the 
flow-Through module. The following equation was utilized to calculate the heat 
transferred to the coolant; 




pQcpfsTj.ff 

60 


(13) 


where 
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Density (g/ml) 

Flow rate (ml / min) 

Specific Heat (J / g “C) 

Average Temperature Differential (“ C) 

c. Average Module Surface Temperature 

The temperature on the surface of the module was measured at 15 different 
locations for each set of experiments. (Side- A and Side- B). The average of these fifteen 
temperatures were calculated for data analyzing purposes as follows: 

\T{i) 

ST^rf=‘^^ (14) 

3. Heat Loss Through the Insulation 

Some amount of heat was conducted through the insulation in either direction 
depending on the power setting and flow rate. Average module surface temperature was 
utilized in calculating the heat losses through the insulation as follows; 

qio.s = ^[sTsurf- T{2)\ Watts (15) 



k = 0.04 W/m “ C Thermal Conductivity of the Insulation 

A = 5.66 in * 6 in = 33.96 in^ = 0.0219 m^ Surface Area of the Insulation 

L = 0.75 in = 0.0191 m Insulation Thickness 

finally, the difference in energy balance is calculated as follows: 

Aq = Pi- q fluid - qioss Watts (16) 
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IV. EXPERIMENTAL RESULTS 


A. HEAT TRANSFER DATA 

A program written in QBASIC code was utilized to obtain experimental data, 
(voltages and temperatures) through the data acquisition ststem. Voltage and temperature 
values were used to calculate the following information in the Tables 6 through 41; 

• Heat suplied to the FTM by each heater, in watts. 

• Total power supplied to the module, in watts. 

• Average module surface temperature, in "C. 

• Average temperatures of the coolant at the inlet and outlet of the module and 
the differential of these two temperature values, in ‘’C. 

• Thermophysical properties, density and specific heat of the coolant, Brayco 
Micronic 889, in g / ml and in J / g “C respectively. 

• Heat dissipated by the coolant, in watts. 

• Heat losses through the insulation, in watts. 

Thermocouples 5 through 19 were used to measure the surface temperatures of 
the module at 15 different locations. These temperature values are plotted in Figures 19 
through 27 and 33 through 37 for both sets of experiments conducted for Side- A and 
Side- B, respectively. 

The graphs of average temperature differential of the coolant between inlet and 
outlet vs. flow rate are shown in Figures-28 and 38 for Side -A and Side -B respectively. 
Logarithmic scales of the same graphs are shown in Figures 29 and 39. Average Module 
surface temperature vs. flow rate is graphed for both Side -A and Side -B in Figures 30 
and 40 respectively. Also logarithmic scales of these graphs are plotted in Figures 31 and 
41.Total power supplied to the module was compared with the sum of heat dissipated by 
the coolant and the heat lost through the insulation and the graphs are shown in Figures 
32 and 42 for both sides of the FTM. 
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Figure 18. Temperature Distributic 
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Table 8. Power and Temperature Data for Side-A at Power Setting 3 
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FTM SIDE-A 
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FTM SIDE-A 

POWER SETTING tf 6 I 70 WATTSI 


Flow Rate 

Heater -1 

350 

POWER 

7.888’ 

450 ! 

POWER 

550 

POWER 

650 

POWER 

750 

POWER 

7.883 ' 

850 ! 

POWER 

7.884 

950 

POWER 

7,875 

7.88 

7.881 

7.884 i 

Heater -2 

7.847 

7.84 

7.842 

7.843 ' 7.844 

7.845 , 

7.837 

Heater -3 

8.787 

7.78 

7.781 

7.783 7 784 

' 7.785 ' 

7.777 

Healer -4 

7,572 

7 564 

7.569' 

7.568 ' 7.573 ' 

7.572 

7.566 

Healer -5 

7 766 

7.759 

7.761 

7.762 ' 7.764 7 

7.764 ■ 

7.756 

Heater -6 

7.828 

7,821 

7.823 

7.824! 7.826 7 

7.826 

7,819 i 

Heater -7 

7769 

7.762 

7763 

7.764 ^76^5 

7.766 7 

7.758 

Heater -9 

7.879 

1%11 

7.8'74' 

rim' 7.877' 

’ 7.878 ' 

77871 

TOT. POW. 

71.109 

70.044 

70.062 

70.075 70.087 

70.0917 

70.022 

Av. Sur. T. 

27.75 

26.49 

25,72 

25.22 ' 24.84 ' 

24.56 

" 24.28 

Tin 

21.12 

20.88 

20.84 

20.83 20.8 

20.81 ' 

20.74 

T out 

28.31 

26.54 

25.5 

24.83 7 24.3 ' 

23.91 ' 

’23.55 

DENSITY 

SP^EAT 

Q fluid 

0.7904 

2,0737 

68.74 

0.7913 
2.0693 
69.51 ' 

0.7917 

2.0669’ 

69.9 

0.792 7 0.7922 ’ 

2.0653 ’ 2.0641 ' 

”70.88' 71.54 7 

0.7924 7 

0.7926 

71 8 ’ 

72.73 

Qins^ 

0.13 

r 0.07' 

0.04 

0.02 0 01 

o’ 

-0 01 ' 

Q total 

68.87 

7 69.58 

69.94 

70.9 71.55^ 

71.8 ^ 

72.72 

DELTAQ 

-2.239 

i -0.464 

-0.122 

0.825 ' 1 463 

1.709^ 

2,698 

DT ins 

" 2.78 

1,47 

1 ' 0.92' 

0.49 ' 0.25 ■ 

0 1 

-0.14 

DT fluid 

7.16 

5.63 

4.72 ' 

'47 3.54 

7 3.09 7 

7 2.75 

T/C-0 

24.97 

25.02 

24 8 

24.73 1 24.59 

' 24.46 

' 24,42 

T/C-1 

’27.04 

26.07 

25,43 

24.99 24,66 

24,4 

24.19 

' T/C-2 

19.93 

19.85 

’ ’ To 97 

20.05 ' 20 

7 20.03 

i’ 20 

T/C-3 

21.19 

'20,99 

' 20.85' 

' 20.88 

20.83 

20.88 

1 20.82 

T/C-4 

28.35 

26.62 

25.57 

24.88 

24 37 

23.97 

23.57 

' T/C-S 

22.96 

’22.48 

22715 

21.97 

21.86 

' 21.8 

I 21.66' 

T/C -6 

24.6 

I 23.9 

23.47 

23.19 

22.99 

22.81’ 

22.68 

' T/C-7 

2609 

25.28 

24.72' 

24.39 

'"’24.11 

23.98 

i 23.72 ■ 

T/C -8 

26.31 

2131 

24,7" 

24.29 

24.04 

23,83 

' 23.57 

T/C-9 

’’ 27.35 

2633 

25.72 

25.34 

25 

24.79 

i 24.57 

T/C-10 

’ 26.61 

25,48 

24.72 

' 24.25 

23.9 

^ 23.66 

23 41 

T/C-11 

2864 

27.38 

26.65 

26.14 

25.71 

25 42 

'25 14 

T/C-12 

28.28 

26 85 

26.03 

25.45 

25 02' 

24,74 

24.44 ' 

T/C-13 

27.91 

26,49 

25.61 

25.04' 

24.61 

24.28 

7 " 23.95 

T/C -14 

27,91 

26.42 

, 25.51 

24.93 

24.52 

24.17 

' 23.85 

T/C-15 

27.48 

26.12 

^ 25.35 

24.78 

24.37 

24.07 

' ' 23,8 ’ 

T/C-16 

29.34 

27.84 

26.92 

26 35 

25.9 

25.55 

25.22 

T/C-17 

30.62 

28.83 

^ 27.69 

27.01 

26.44 

26.03 

'2’5.65 

T/C -18 

31.26 

29.61 

1 ’ 28.53 

27 81 

27.25 

26.86 

26.45 

T/C -19 

30.83 ■ 29.07 

i ' 28706 

27 34 ■ 26.82 

26.42 

26.04 


Table 11. Power and Temperature Data for Side-A at Power Settmg 6. 
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FTM SIDE-A 



POWF.R 

SETTING 

#7/80 WATTS! 



Flow Rate 

350 ' 

POWER 

450 

POWER 

550 

650 

750 

850 

950 ‘ 

POWER 

9.001 

POWER 

9.002 

POWER 

POWER 

POWER [ 

9.007 ' 

8.992 

9 001 

9.001 

9,002 1 



8 946 

8.956 

8.955 

8.956 ; 

8.957] 

8.957 


8.892^ 

8.877 

8.88*7 

8.887 * 

8.888 

8.888 1 

8.888 


8.659 

8.648 

8*65*2 

8.65*4** 

8.65* 

8.644' 

8.645 


8.868 

8 855 

8.864 

8.863 ■ 

8.864 1 

8.864 J 

8.864 




' 8.934 

8.934* 

8 935 1 

8 935 1 

8.936' 

Heater -7 

8.869 

8.859 

8 864 

8.864*' 

8 864 1 

8.865 

8.865 

1 Heater -8 

8.874 

8.866 

8.87 

8 87 

8 ^ 

8.871 1 

* 8.872^ 

'lleater -9 

8.996 

8,987 

8.992 

8.991 

8.993 

8 993 

8.993 ; 

TOT. POW. 
Av.Sur.f. 

80.064 

28.85 

79.954 

26.6 

80.021 

26.52 

80.019 

25 99 

80.021 ’ 

25.54 

80.019 1 

25 15] 

80.021 

24,9 

Tin 

21.31 

20.99 

20.97 

20.97 

20.85 ' 

20.92' 

20.83 

T out 

29.45 

27.36 

26 23 

25.54 

24*88 

24.4’ 

2401 

DENSITY 

0.7899 

0.7909 

0J914 

0.7916 

0 792 

0.7921 1 

0.7923 

SP.HEAT 

2.0767J 

2.0713 

2.0688 

2.0672 

2.0655 ■ 

2.0646' 

2.06357 

Q fluid 

77.89 1 

78.2*6' 

78.94 

81.02 i 

82.41 ' 

80.62 ‘ 

82.‘32l 

Qins ; 

0.14J 

6.05 

0.05 

0.03 i 

*0.01* 

0 ' 

' -0.01 

Q total 

78.03 : 

7831 

78.99 

81.05 

82.427 

80.62 

8231 ! 

DELTAQ 

-2.034 

-1.644 

“ -1.031 


; 2.399 

0.601 

r * 2.289 

DTins 

3.14 

*1.19 

.1 .02' 

0.6 

' 0.24 i 

-0.09 

1 -0.16 

' DT fluid ~ 

8.14 

5.31 

5.25 

4.49 

; ■ 3.97 

3.47 

i 

THFRMOrODPI 

K RFAD' 

INC.S 

T/C-O 

25.71 

25.41 

25.5 

25.39 

I-2 ^ 

25.24 

1 25.06 

T/C-1 

28.02 

26.22 

26.16 

:* * 25'7' 

I 25.31 

25.05 

1 * 24.8' 

T/C-2 

19.99 

19.96 

20.04 

20.06 

l_ 19.97 

20.12 

j 19.93 

T/C-3 

21.4 

21.13 

21.03 

21.03 

1 21.04 

1 20.94 

20.91* 

T/C-4 

29.54 

26.44 

26.28 

25.52 

25.01 

2441 

24.09' 

T/C -S 

23.45 

22.53 

22.56 

22.25 

22.13 

' 2199 

]_ 21.86 

T/C-6 

25.26 

2402 

■ 23794 

23.63 

23.41 

23.22 

*23*.04 ! 

T/C-7 

27.02 

25.45 

' 25.41 

25.07 

24.76 

24.48 

i *24.23* ’ 

T/C-8 

27 18 

25.44 

25.39 

24.92 

24.62 

* 24.3 

■ *24.05 ' 

T/C-9 

28.4 

26.57 

26.58 

26.08 

2576 

25.4 

^ 25.27 

T/C-10 

27.51 

i ' 25.47 

25.35 

24.86 

' 24.46 

24.12 

i 23.91 i 

T/c-n 

i 29.88 

i 27.59 

27.55 

26.96 

' 26.53 

26.2 

I 

T/C-12 

i 29.42 

1 26.94 

26.82 

26.3 

j* 25.76 

25.28 

I 2604 ! 

T/C-13 

I" 29.01 

I * 26.46 

■ 26.36 

25 79 

1 25.2*3 

; * 24.86 

i 24.54 1 

T/C-14 

1 29 

26.38 

26.31 

■ 25.71 

25.18 

' 24.72 

24.43 

T/C -is 

L 28.53 

26.19 

26 

25 51 

25,04 

24*6* 

i 24.35 

T/C -16 

30.69 

28.01 

27 93 

27.32 

26.76 

; 26.28 

. 26 

T/C-17 

32.08 

28.89 

28 76 

28.03 

^ 27.34 

' * 2683 

26,44'’ 

T/C-18 

32.94 

29.82 

29.67 

2891 

28.32 

27.79 

27 44 

' T/C -19 

32.34 

29.23 

29.14 

28.44 

27.76 

* 27.2*5* 

, 26.95 


Table 12. Power and Temperature Data for Side-A at Power Setting 7 
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Figure 25. Temperature Distribution on the Module for Side -A, at Power Setting 8 
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Figure 26. Temperature Distribution on the Module for Side - A, at Power Setting 9. 
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COOLANT TEMPERATURE DIFFERENCE vs FLOW RATE I 
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LOG-LOG PLOT OF TEMPERATURE DIFFERENCE vs^ FLOW RATEI 
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I AVERAGE SURFACE TEMPERATURE vs FLOW RATE 












TOTAL POWER vs HEAT RATE COMPARISON! 













FTM SIDE - B 


POWER SETTING i 

1 1 (20 WATTS)- 

Flow Rate 

350 

450 550 

650 

750 

850 

950 


POWER 

POWER POWER 

POWER 

POWER 

POWER 

POWER 


2.249 

2.25 2.247 

2.247 

2.247 

2.248 

2.248 



2.239 : 2.236 

2.236 

2.236 

2 236 

2.237 

Heater -3 

2.221 

2.222 2.219 

2.22 

2.22 

'2.22 

2.22 

Heater -4 

2 185 

2.186 2.184 

2.184 

2.184 

2184 

2.185 


2.215 

2.216 2.213 

2.213 

2.213 

2.214 

2.214 



2.234 2.232 

2.232 

2.232 

2.232 

2,233 


2.218 

2.219 2.216 

2.216 

2.217 

2.217 

2.217' 

Heater -8 

2.218 

2.219 2.217 

2.217 

' 2.217 

2.218 

2.218 

Heater -9 

_ 2.249 

2.25 ^_ 2.247 

2.248 

2.248 

2.248 

2.248 


22.46 

^2^5 i Tt75^ 

21.54 

21.36 

'21.^ 

21.12 

^Tin “ 

20.92 

20.81 : 20.7 

20.59 

20.52 

20.5 

20.41 

Tout 

22,92 

22,34 2'l.9S 

'21.65 

21'45 

" 21.32' 

'21.17' 

DENSITY 

0.7928 

0.793 0.7932 

0.7934 

0 7935 

0.7936 

0,793^7 

QOuid 

19,07 

l'8.74 i 18^1* 

“ 18^5 

1898 

18 96 

19.64 

Qina 

-0.06 

-0 08 i -0 08 

-0.09 

-0 09 

-0,1 

-0.1 

0 total 

19.01 

18.66 18.63 

18.66 

18.89 

18.86 

19.M 

DELTA Q 

-1.016 

-1.375’ -l'.'3"81 

-1.353 

-1.124' 

-1.157 

-0.48 

DTim 

-1.34 

-166’ -1.78 

-1,9 

-1 98 

-2.12 

-2.15 

DT fluid 

1 1 94 

1 51 ' 1.26 

1.04 

1 

0 85 

0.77 

THERMOCOUPl 

J:READINGS 

T/C-0 

1 19.23 

19.4 i 19.47 

19.27 

19 44 

1953 

19.5 

T/C-1 

' 2525' 

24 88 24.68 

24.43 

■' 24.27 

24.18 

24,07 

T/C-2 

23,8 

23.71 23.53 

23.44 

23.34 

23.34 

23.27 

T/C-3 

21 02 

20 89 ' 20.77 

20,74 

20.58 

20.54 

20.42 

T/C-4 

22.96 

22.4 ! 22.03 

21.78 

2158 

21.39 

21.19 

T/C-5 

21 99 

21.64; 21.4 

21.23 

21.05 

20.98 

20.88 

T/C-6 

21.2 

21 ' 20 84 

20.79 

' 20.7 

20.61 

20.5 ! 

T/C-1 

21,92 

21.64; 21.'46 

21.32 

^ 2143 

21 

20,93 1 

T/C-8 

22.14 

21.82 21.59 

21.42 

1 21,24 

21.12 

21.14 

T/C-9 

22.07 

21.77 2'r.49 

21.31 

i' 21.13 

20.99 

20.91 1 

T/C-IO 

21.87 

21.58 ■ 21.34 

21.17 

21 

20,83 

20.8 

T/C-11 

22.1 

21.71 ^ 21.44 

21.29 

21.12 

20.95 

20.92 j 

T/C -12 

22.71 

22.26' 21.93 

21 72 

21.52 

21.36 

21.22 

T/C-I3 

22.94 

22,42 ’ 22.09 

' 21.85 

21,65 

21.49 

21.37 1 

T/C -14 

22.8 

22.3 21.94 

21.72 

21.58' 

21.44 

21.29 

T/C -IS 

22 54 

'22.12,_ 218 

21.56 

2l'35 

21.24 

21.05 

T/C-16 

22,61 

22.16' 21.82 

21.56 

21.39 

21.26 

21.11 

T/C-17 

23.28 

22.72 " 22.35 

22.03 

2T82 

21.64 

21.52' 

T/C-18 

23 5 

22,94 ■ 22.5 

_22.24 

2199 

21.82 

21-69 1 

T/C -19 

23.18 

; 22 65 , 22.23 

1 21.94 

; 21.72 

21.55 

21.44 


Table 15. Power and Temperature Data for Side - B at Power Setting 1 
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; of Coolant Inlet and Outlet Temper; 
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B. FORMULATION OF THE CORRELATIONS 

The correlation calculations are based on some assumptions made for the free 
flow area, Af and the log mean temperature difference, ATlm. The assumption for the 
free flow area is made due to the uncertainty of the iimer dimensions of the FTM. This 
assumption is made for the width of the 3 parallel passages in order to calculate the 
number of channels in a passage shown in Figure 7 and detailed in Figure 42. 

Heating of the FTM would not fit either the "constant surface heat flux" case, or 
the "constant wall temperature" case. Therefore the log-mean temperature difference was 
defined for the experiments as follows: 



(7’(i9)-r) -(7'(5)-r,„) 


The rest of the necessary data is available in order to be able to calculate the 
Reynolds and the Stanton numbers. The equations utilized for the calculations of these 
quantities are as follows: 

(18) 

= - (19) 

pQcpATiM 

The hydraulic diameter is calculated by using the following equation and the 
dimensions defined for the fins, in Figure 6 and drawn in detail in Figure 42. 
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( 20 ) 


D, = ^ = 4 X 0 046 X 0.06 ^ 0.05208/w = 1.323 x IQ-^m 

P 2 X (0.046+ 0.06) 

The free flow area of the coolant through the three parallel passages of the FTM 
is evaluated with respect to the assumption made for the passage width ( Figure 42) 
which is assumed to be 1.75 in. Since the fin density is known the free flow area is 
calculated as follows; 

Aj = bxn (21) 

= 0.06 X 0.046 X 34 = 0.09384^2 = 6.0542 x lO-^m^ 


where, 

h: Height of a Channel. 

h'. Width of a Channel. 

n: Number of Channels in a Passage in the Transverse Direction. 

The average velocity of the coolant through the channels is calculated by using 
the following equation: 




Q 

Af 


( 22 ) 


Finally, the necessary information to evaluate the Reynolds and the Stanton 
numbers are available and the results are tabulated in Table 22 and plotted in Figures 43 
through 58. 

The relationship between Reynolds and Stanton numbers is formulated by 
utilizing the logarithmic plots. The results for each power setting are tabulated in Tables 
20 and 21, for both sides of the FTM. 
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Figure 42. 







! CORRELATION (St vs. Re) 

TOTAL POWER 
(Watts) 

St = 0.100*Re-®“‘ 

30 

St = 0.108*Re-"™ 

50 

St = 0.111*Re^™ 

60 

St = 0.110*Re-<'’*^ 

70 

St = 0.120*Re-“*”’ 

80 

St = 0.123*Re-®“ 

90 

St = 0.130*Re-®*^’ 

100 


Table 20. Correlations for Side - A. 


1 CORRELATION (S/vs./?e) 

TOTAL POWER 
(Watts) 

St = 0.084*Re-'’‘'“ 

20 

St = 0.091*Re-“^'‘' 

40 

St = 0.098*Re^^'’ 

60 

St = 0.084*Re-"'*^ 

80 

1 St^O.lOPRe'"^"’ 

100 


Table 21. Correlations for Side - B. 
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SIDE-A 

30 watts 


Flow Rate 

350 

450 

550 

650 

750 

850 

_ 950 

I.ogmeanT 

0938 

0 056 
2.03 

0.901 

0.0454 

2.58 

0.867 

0837 

0.796 

0.804 

0.769 


317' 

3 66 

4.19 

4173' 

5.29 

Kinm.Visc7 

0.096 

0.000063 

0.124 

0.000064 

0.151 

0.000064 

0.179 

d.00TO65 

0.206 
0.000065 i 

0.000065“ 

0,262 

0 000066 




70 watts 





Flow Rate 

350 

450 

550 

650 

750'; 

850 

950" 

LogmeanT 

2.162 

203 

1.866 

. 17361 

1.686 ‘ 

1.634' 

1.577 

St# 

0.0573 

0 0468 

0.0417 

0.038 ' 

0.0339 

0.0309 ■ 

0,0286 

Re# 

2 36 

2 87 

3.4 

3,95 j 

4 48 ■ 

' 5 04" 

_5.58 

llmMn 

Kinm.Visc. 

0.096 

0.000054 

0.124 

0000057 

0.151 
n 000059 

0 r79 I 
o.ooo66 

0.206 ' 
0000061 ‘ 

0.234' 
0 000061 ' 

0 000062 




100 watts 





Flow Rate 

350 

450 

550l 

~ 650 

750 

850 i 

950 

LogmeanT 

2.996 

2.723 

2 551 1 

;_^4^7 ^ 

2.375 ^ 

2.283 

2.232 

Re'' 

2.54 

3,03 

1.751 

r °'°4r 

4^62 ■ 

^,n 

5.7 


0.096 

0,124 

I 0.151 

0.I79' 

0.206 ' 

0.234 

0.262 

Kinm.Visc. 

1 0.00005 I 

0.000054 

; 0.000056" 

0 000058 ' 

0.000059 ' 

0 00006 

0.000061 




S [DE- 

B 




60 Watts 

Flow Rate 

350 

450 

550 

650 

750 

850 

950 

LogmeanT 

1 416 

1.194 

1.057 

0.95 

0,848 ' 

; 0.779 

0.71 

Re# 

2.29 

2.81 

3.32 

3 88 

4.39 

4.94 

5.48 

L memi 
Kinm.Visc. 

0.096 
' 0.(X)0055 

0 000058 

0.151 

000006 

0,179 

' 0.000061 

0.206 

0,000062 

0.234 

0.000063 

0,262' 
0 0000631 




80 Watts 





Flow Rate 

350 

450 

550 

650 

750 

' 850 

950 

LogmeanT 

1.913 

0.0729 

2.4 

1.545 

0.0703 

2.91 

1.265 
0.0703 

1 135 
0.0663 
3.98 

1 1.064 

1 0.0613 

U 4.5 

1 0.9 

' 0.064 

5.05 

; 0.808 
[ 0.06381 

5.^ 

Re# 

Kinm.Visc. 

0.096 
0 000053 

0.124 

0,OOOU56 

0.151 
0.000058 

0.179 
’ 0.00006 

0.206 

0.000061 

0.234 

0.000061" 

0.2621 

0.000062 




100 Watt! 

i 




r Flow Rate 

350 

450 

550 

650 

?«)" 

850 


LogmeanT 

_ 2 281 

1 842 

1 641 

’ 1 456 

i 1.332 

; 1 146 

L — 

Re# 

25 

2.9S 

"" 7.49 

i ^ T05 

~" 4TS4 

5.09' 

T 5.62 

U mean 

’ 0 096 

‘ 0.124_ 

' _ 0.151 

0.179 

0.206 

0.234 

0.262 

Kinm.Visc. 

0 000051 

0.000055 

0.000057 

1 0.000059 

0.00006 

0.000061 

0.000062 1 


Table 22, Reynolds and Stanton Number Data for Both Sides of the FTM 
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STANTON # vs REYNOLDS #| 


SIDE - A (30 W 



Figure 43. Linear Plot of Stanton vs. Reynolds Number for Side -A, 
at Power Setting 2. 



Figure 44. Logarithmic Plot of Stanton vs. Reynolds Number for Side -A, 
at Power Setting 2. 







STANTON # vs REYNOLDS # 

SIDE - A (70 W) 



Figure 45. Linear Plot of Stanton vs. Reynolds Number for Side -A, 
at Power Setting 6. 


I STANTON # vs REYNOLDS # 

SIDE - A (70 W) 



001 


REYNOLDS NUMBER 


Figure 46. Logarithmic Plot of Stanton vs. Reynolds Number for Side -A, 
at Power Setting 6. 


10 
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Figure 49. Linear Plot of Stanton vs. Reynolds Nuinber for Side - A, at Power 
Settings 2,6 and 9. 



Figure 50. Logarithmic Plot of Stanton vs. Reynolds Number for Side - A, at I 
Settings 2,6 and 9. 








at Power Setting 3, 



Figure 52. Logarithmic Plot of Stanton vs. Reynolds Number for Side -B, 
at Power Setting 3. 
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STANTON # vs REYNOLDS # 

SIDE - B (80 W) 


0.074 
0.072 i- 



0.064 


0 062 


Figure 53. Linear Plot of Stanton 
at Power Setting 4. 


REYNOLDS NUMBER 
vs. Reynolds Number for Side -B, 



at Power Setting 4. 
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Figure 57. Linear Plot of Stanton vs. Reynolds Number for Side - B, at Power 
Settings 3,4 and 5. 



figure 58. Logarithmic Plot of Stanton vs. Reynolds Number for Side - A, at 1 
Settings 3,4 and 5. 
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V. DISCUSSION AND CONCLUSIONS 


Both sets of experiments conducted for Side-A and Side-B clearly showed the 
effectiveness and high cooling capacity of the FTM. 

A. THE EXPERIMENTS CONDUCTED FOR SIDE-A 

The thermocouple readings plotted in Figures -18 through 26 show the general 
trend of temperature rise along the first and third rows and temperature decrease in the 
second row of the flow path. 

The decrease in the surface temperatures of the module, as shown in Figures -18 
through 26, was observed more clearly during the changes of flow rates from 350ml/min 
to 450 ml/min and from 450 ml/min to 550 ml/min at all power settings. For the rest of 
the flow settings the decrease in surface temperatures were almost equal and actually the 
decrease varied logarithmically. The lowest temperature was measured at thermocouple 
5, which was the first reading taken from the surface of the module. This low temperature 
occurred due to the which fins existed on the path from the inlet to that thermocouple 
location. On the other hand, the temperature on the surface of the module was measured 
at thermocouple location 18, which was before the last thermocouple affixed just after 
heater-F (Figure -12). The slope of temperature rise across the heaters is almost equal for 
the heaters A,B,C and E. This slope was steeper than the others, that is temperature rises 
across these heaters were more significant than heater D and F. 

On the other hand, the temperature was observed to decrease across heater-D for 
all of the power and flow settings. This was explained by the counter-flow heat exchanger 
effect occurring between the inlet and the bend, between the second to the third rows. 
Significant negative slope was observed between thermocouple readings 11 through 15 
due to this counter-flow heat exchanger effect. 
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The temperature difference between inlet and outlet of the module, showed a 
logarithmic decrease as the flow rate increased (Figure -27). A linear relationship was 
observed when the same graph was plotted on log-log scale (Figure -28 a and 28 b). 
These linear lines were observed to be parallel which helped in determining the 
correlation in terms of Stanton number and Reynolds number. The average surface 
temperature of the module showed the same trend of logarithmic decrease as the flow rate 
increased. 

The heat dissipated by the coolant for all power settings was observed to be 
increasing as the flow rate increased (Figure 31). A tendency of higher heat dissipation 
level than the total power supplied to the heaters was observed. This arose in 33 of 63 
experiments ( 50%) conducted for Side-A, with the highest increase of 10% (Power 
setting 2). Similar results were generated in the previous study and in the experiments 
conducted by the Crane Division of Naval Surface Warfare Center. This occurred in 60 of 
80 runs (75%), during the previous study and 16 of 23 runs (69.6%) during the 
experiments conducted by NWSC [1,2]. According to the guidelines presented in the 
previous study related to the coolant specific heat, Delsen Laboratories data was used in 
calculations, in order to be able to acquire better results than before. However as 
mentioned earlier in 33 of 63 runs conducted for Side-A the same results as before were 
observed. 2.3% of the discrepancy in the heat dissipation level was explained by the 
uncertainty , calculated for q , in Appendix - B. However there still some amount of 
energy difference existed between the power into heaters and the heat dissipated by the 
coolant. 

The general trend was observed as follows: The heat dissipation level was less 
than the total power supplied into the heaters at lower flow rates and it kept increasing as 
the flow rates increased. Thus the amount of heat dissipated by the coolant was more than 
the total power supplied by the heaters at high flow rates. 
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The criteria defined for steady state achievement, was assumed to be the same for 
higher flow rates. However some amoimt of fluctuations were observed at the inlet 
temperature of the coolant as the flow rate increased. This resulted in noticeable 
difference between the temperature differentials, even though steady state was assumed to 
be achieved. This might be the cause of higher heat dissipation levels observed at higher 
flow rates. 

The general form of the correlations obtained by formulating the relationship 
between the Stanton and the Reynolds numbers was as follows: 

In (St) = m* In (Re) + n 


The coefficient A, increases as the power increases and the maximum difference 
between the values of exponent B, is about 0.08 which is assumed to be in a reasonable 
range (Table 20). The increase in the coefficient A, can be explained by the Prandtl 
number effect of the coolant. 

B. EXPERIMENTS CONDUCTED FOR SIDE-B 

The decrease in the surface temperatures of the module was observed more clearly 
during the changes of the flow settings at lower flow rates as observed on Side-A. 
Cooling capability of the FTM was more significant during the first two flow settings 
which was the same result obtained from the first set of experiments. The same 
logarithmic trend, in the decrease of module surface temperatures was observed as the 
flow rate increased. The lowest and the highest temperatures consistently occurred at 
thermocouples 6 and 18 respectively. A sharp decrease in temperature was observed 
between thermocouples 5 and 6 which could be explained by the existence of the fins 
along this path. 
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The same trend of logarithmic decrease was observed in the average surface 
temperature vs. flow rate graph. Unlike Side-A, the heat dissipated by the coolant was 
always less than the total power supplied by the heaters. 

The relationship between Reynolds and Stanton numbers was formulated in the 
same way as for Side - A. The coefficient A was also increasing as the power rate 
increased, which verified the Prandtl number effect of the coolant. The maximum 
difference in the values of the exponential B, is about 0.18. 
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VI. RECOMMENDATIONS 


In order to be able to completely interpret the data obtained throughout this and 
the previous studies [1,2] knowledge of the internal configuration of the FTM would be 
required. This would also provide a valuable information in acquiring the correlation by 
making the corrections to the ffee-flow area. Unfortunately due to the proprietary nature 
of the FTM, it was not possible to cut it open. 

The experiments can be repeated using water as the coolant of the system to 
compare the data and to determine the Prandtl number dependence on the effectiveness of 
the module. 

The correlations in terms of Colburn j factor can be calculated using the data and 
the results can be compared with the ones obtained in a previous study [14]. 

The surface temperature patterns at various power and flow settings can be 
visualized by using thermochromic liquid crystals. The temperature data obtained in this 
and in the previous study can be compared with the pictures of these images [1], 
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APPENDIX A. 


SAMPLE CALCULATIONS 


Sample calculations are based on power setting # 9 (100 Watts) of the first set of 
experiments for the flow rate of 350 mVmin. The results are tabulated in Table 13 and the 
procedure can be summarized as follows; 

The current passing through each heater is calculated by using Equation 8 as: 

where i is the index for heaters and precision resistors from 0 to 8. Thus for 
power setting #9 current passing through heater # 8 is calculated as; 

U„(8) 

«S).W8) = ^ 

and using Equation 9, the power supplied to the heater # 8 is calculated as: 

P/.(8) = Va(8) X 4(8) = 15.89 ♦ 0.707 = 12.23 Watts. (9) 

from Equation 11 heat flux generated by each heater is evaluated as: 

q" = = Watts/cm^ (11) 

where A = 1.94 cm^ 


average of the differential between coolant inlet and outlet temperatures is 
calculated by using Equation 12 as follows: 
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Tdiff-'Z- 


( 12 ) 


The thermophysical properties of the coolant, that is density and specific heat, are 
calculated by using Equation 1 and 3 respectively as follows: 


p =-8.8216 xlO-’x 7^,^-7.9251 xr„v^+0.81055 = 0.7894 g/ml. (1) 


Cp = (3.528 x\0^xTl,g + \.047 x lO'^ x Tavg + 0.4695) X 4.184 


= 2.0791 


J/g°c 

where Tmg 


Ti(.4)-U3) 

2 


= 25.92'’C 


(3) 

(2) 


The volumetric flow rate, which is chosen to be 350 ml/min is calculated by using 
the following fourth order polynomial (Equation 5) obtained from the flowmeter 
calibrations: 


2 = 2.7178xl0-5a''-6.02xl0-3a3+0.4224a^+4.1837a-8.073 ml/min (5) 

where is the flowmeter reading in percent, %. Therefore heat dissipated by the 
coolant can be calculated by using Equation 13 as follows: 

q/iuid = = 97.75Waffs (13) 

The average module surface temperature, which is utilized in heat loss calculation, 
is calculated by Equation 14 as follows: 

Tsu^ = ‘^^ = 30.l6‘’C (14) 
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The heat loss through the insulation is calculated by using Equation 15 as follows: 


qioss = ^[T.urf-T{2)\ 


(15) 


= 16 - 24.27) = QAlWatts 

finally the difference in energy balance is calculated as; 

Aq = P,- qjiu,d - qion = 100.025 - 97.75 - 0.47 = 1.805M'arrs ( 16) 


The Reynolds and Stanton numbers are calculated by using Equations 18 and 19 
respectively. Mean velocity of the fluid through the channels for the flow setting 350 
ml/min is evaluated as follows: 


f/m 


Q 

Af 


(22) 


ml s/ \m^ w 1 min.. 


min-^M0«m/ 60 sec6.054 X 10-5 


= 9.63cmls 


The log mean temperature difference, used in Stanton number calculations is 
calculated for power setting # 9 as follows: 


[r(i9)-r„„,]-[r(5)-7’,„j 

ln[[T(19)-7-„„,)^[7-(5)-r,„]] 


( 17 ) 


(34.4 - 31.03) - (23.47 - 20.82) 

In [(34.4 - 31.03) - (23.47 - 20.82)] 


= 2.996''C 
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Using Equation 18, the Reynolds number is calculated as follows: 


0.096x1.323x10-3 

0.00005 


(18) 


and finally, the Stanton number is evaluated by using Equation 19 as follows: 


" pQcpATiM 0.7894 x 2.0791 x 350 x 2.996 


(19) 



APPENDIX B. UNCERTABVTY ANALYSIS 


There was some uncertainty in measurements taken during the experiments. These 
uncertainty values were either evaluated or estimated for all of the variables mentioned in 
sample calculations and summarized in Table 23. 

The uncertainty of a function in the form, / (x^ ,xj = x, * x^ is evaluated by 
utilizing the following equation: 




where; 

X, and Xj 
8 / 

5x, 

8x: 


experimental variables. 

Total Uncertainty in Function f. 
Uncertainty in Variable x,. 
Uncertainty in Variable x^. 


(23) 


The uncertainty calculations were performed for power setting # 5 and flow setting 
of 950 ml/min. Where: 


p = 0.7929 g/ml 

V^ =12.288 Volt 

Vp= 0.275 Volt 

v = 6.3 X 10-’/wV5 

U„ = 0.262 m/s 

Re = 5.48 

Cp =2.0607 J/(g°C) 

1.332‘’C 

St = 0.0291 

fj = 20.56“C 

7’4 = 23.00°C 


7.,^ = 21.78‘■C 

Rp = 0.5009 Ohm 


Therefore the uncertainty in 

power supplied to the heaters can be calculated using 

Equations 9 and 7: 



P, = V,xh 


(9) 

Vn 


(7) 
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1 VARIABLE 

UNIT 

UNCERTAINTY 

REFERENCE I 

1 Rp 

Ohm 

0.0001 

Calibration of Precision Resistors | 

V„&V, 

Volt 

0.008 %+0.0007 

Integrating Voltmeter Manual 1 1| 

Ip 

Amp. 

0.00263 % 

Calculated Eqn.(24) || 

P, 

Watt 

0.00263 % 

Calculated Eqn.(25) 

T 

°C 

0.12‘’C 

Thermocouple Calibration 

Tdiff 

”C 

0.078% 

Calculated Eqn.(26) 

P 

g/ml 

0.005 

Estimated 

Cp 

J/g°C 

0.005 

Estimated 

|v 

mVs 

10® 

Estimated 

a 

% 

1.5 

Experimental Observations 

Q 

ml/min 

2.1% 

Interpolation 

A 


10 %(0.00219 m"') 

Estimated 

L 

m 

0.003m 

Ruler Resolution 

T surf 

"C 

0.078 % 

Calculated Eqn.(26) 

'ATlm 

°C 

1.52% 

Calculated Eqn.(36) 

h 

m 

10® 

Estimated 

b 

m 

10® 

Estimated 

n 


2.94 % 

Estimated 

K 


1.08% 

Calculated Eqn.(30) 

Ap 


3.13% 

Calculated Eqn.(33) 

D, 

m 

1.53% 

Calculated Eqn.(29) 


m/s 

3.77% 

Calculated Eqn.(35) 

qi.ss 

Watt 

23.3 % 

Calculated Eqn.(28) 

q 

Watt 

2.34 % 

Calculated Eqn.(27) 

Re 


4.37 % 

Calculated Eqn.(37) 

St 

-- 

3.6% 

Calculated Eqn.(3 8) 


Table 23. Uncertainty Values of the Variables 
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(24) 


+ =2,63xl0“3 




(25) 


=±0.158»^ 

The uncertainty in average temperature differential from Equation 12: 

Tdiff= T4 - Ti 


1(57:4)2 +(1^)2 =0.78 (26) 

T 4 Tj 

8Td,ff=+0.0\^6^C 


The uncertainty in heat dissipation satisfied by the coolant is evaluated from 
Equation 13 by utilizing the above results as follows: 


|(^)2+(5e 2 ^ ^ 

qfl.id ^^Td^ 

5,^„=±l AW 


= 2.34% 


(27) 


The uncertainty in the flowmeter reading was +1.5% with respect to the 
experimental observations. Therefore the uncertainty for this flow setting (950 ml/min) is 
calculated simply by interpolating according to the flowmeter calibration, and found to be 
2.1%. 
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The uncertainty in heat losses through the insulation is calculated as follows: 


^ = 23,3% (28) 


The uncertainty calculations based on Reynolds and Stanton numbers are evaluated 
by calculating the uncertainties in the following parameters as follows: 

The uncertainty inliydraulic diameter: 


Dh 


(29) 


5Ac 5P 
Ac P 


^(§)^+(y)^ =0.0108 


since A ^ = h*b 
and P = 2*{h + b) 

Therefore the uncertainty in hydraulic diameter is 


^ = 7(0.0108)2+(0.0108)2 = 0.0153 


(30) 

(31) 

(32) 


The uncertainty in free-flow area is evaluated as follows; 

=0,03,3 

where A,^h*b*n 


(33) 

(34) 
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By using the result obtained from Equation (34), the uncertainty in the mean 
velocity of the coolant through the channels is calculated as follows: 


U„ 


Q 


(35) 


= 7(0.021)2 +(0.0313)^ = 0.0377 

where Um = ^ 

Af 


( 22 ) 


Finally, the uncertainty in the log-mean temperature difference is evaluated as 
follows: 


5(ATim) 

+ 2 (^+ 2 (^ 


(36) 


= j2(:^)'+2(^)^ + 2(:^^ 


4o.56^ 


The above results can be utilized to calculate the uncertainties in the Reynolds and 
the Stanton numbers as follows: 


^= 0.0437 




- 9 X 0.02367 = 0.0213 is the uncertainty in total power. 
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APPENDIX C THERMOPHYSICAL PROPERTIES OF BRAY CO MICRONIC 889 


The properties of the coolant, Brayco Micronic 889, are tabulated according to the 
data provided by the manufacturer, Castrol [3], as follows: 


TEST(ASTM) 

DESCRIPTION 

RESULT 

D-287 

Specific Gravity 

0.08 

D-445 

Kinematic Viscosity, (10® mVs) 



@ 100"C 

1.7 


@ 40 ®C 

5.2 


@ - 40 “C 

260 


@ - 54 ® C 

1,200 

D-877 

Dielectric Strength, (KV) 

35 

D-1169 

Resistivity, 25® C (77 ® F), (ohm-cm) 

1.5-10'® 

D-287 

Density, (g / ml) 



@ 0®C (32®F) 

0.811 


@ 20®C (68®F) 

0.794 


@ 40®C (104®F) 

0.777 


@ 100®C (212®F) 

0.723 


@ 160® C (320®F) 

0.661 

D-2766 

Specific Heat ( cal / g ®C ) 



@ -18®C (0®F) 

0.49 


@ 10 ® C (50 ® F) 

0.52 


@ 38 ®C (100®F) 

0.54 


@ 93 ®C (200 ®F) 

0.58 
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APPENDIX D. QBASIC PROGRAM CODE 


This Program was written in order to establish an interface between HP 3 8 52A 
Data Acquisition / Control Unit and the Personal Computer. The program code is in 
QBASIC and it is located under the program group, "Thesis Items" in the PC. After 
operating the HP3852A, the program can be run by pressing F5 or by choosing "Run" 
from "File" menu while the file named "FTM.BAS" is open. Temperature and Voltage 
values acquired by HP 3852A are utilized throughout other calculations. 

CLS 

COLOR 14, 9 

PRINT " Program for National Instruments PC2 IEEE-488 BOARD to HP-3852" 

PRINT " This uses the Universal Language Interface" 

PRINT " ULI.COM must be run prior to running the program" 


e of the program is FTM.BAS" 


'Preparing the interface between program and PC2 board 
OPEN "GPIBO" FOR OUTPUT AS #1 
OPEN "GPIBO" FOR INPUT AS #2 

'Initializing the bus and reset to default parameters 
PRINT #1, "ABORT" 

PRINT #1. "RESET" 


PRINT # 1, "GPIBEOS CR LF" 'Establish end of string character. 

PRINT # 1, "CLEAR 9" 'Clear device 9 which is HP-3852 . 

STARTIT: 

PRINT #1, "REMOTE 9" 'Place it into Remote mode to accept 

commands. 

PRINT #1, "OUTPUT 9; RESET 700" 'Reset voltmeter in slot 7. 

PRINT #1, "OUTPUT 9; USE 700" 'Use voltmeter in slot 7. 

PRINT #1, "OUTPUT 9; CONF TEMPT" 'Configure for temperature, type T. 

************* BEGIN TO TAKE DATA*»**************** 

DIM TEMP(20), V( 18), R(9), P(9) 'Defining dimensions for variables. 
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INPUT "ENTER FLOWRATE IN mL/min >", Flow 

FOR 1 = 0 TO 19 'For/Next Loop for channels; All 20 

channels of 

' the HP 44708A relay multiplexer. 

CHNLS = STR$(I + 300) 

MSGS = "OUTPUT 9;MEAS TEMPT" + CHNLS 

'Command to measure temperature 


on channel I. 


PRINT # 1, MSGS ' Sending MEASURE command to HP-3852. 

PRINT # 1, "ENTER 9" 'Call enter routine on interface 

board. 

INPUT #2, TMPS 

TEMP(I) = VAL(TMPS) 'Data returning into TMPS. 

PRINT "THERMOCOUPLE "; (I);" = "; TEMP(I) 

NEXT I 


PRINT # 1, "OUTPUT 9;CONF DCV" 'Configuring for DC Volt measurements. 

PRINT #1, "OUTPUT 9;USE 700" 


FOR I = 0 TO 17 '18 channels are used for Heaters 

and 

' Precision Resistors. 

CHS = STR$(I + 200) ' Slot number for DCV measurements 

is 2. 

MESSS = "OUTPUT 9;MEAS DCV" + CHS 'Measuring voltages on channels 
200+1. 

PRINT #1, MESSS 

PRINT # 1, "ENTER 9" 'Calling ENTER routine. 

INPUT #2, VOLTS 

V(I) = VAL(VOLTS) 'Data is stored to V(l). 

NEXT I 
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PRINT #1, "OUTPUT 9;BEEP'' 
PRINT #1, "CLEAR 9" 

PRINT #1. "LOCAL 9" 

R(0) = .4876 
R(l) = .4877 
R(2) = .4877 
R(3) = .4878 
R(4) = .4876 
R(5) = .4876 
R(6) = .5008 
R(7) = .5 
R(8) = .5009 
TMPT = 0! 
temperature. 


'Precision Resistors in Ohms. 


'Calculating the module average surface 


DO WHILE I < 20 


1 = 1+1 

TMPT = TMPT + TEMP(I) 

LOOP 

T SURF = TMPT /15 
PRINT "AVG SURF. TEMP"; 
PRINT USING "##.##"; TSURF 
FOR I = 0 TO 8 
P(I) = V(I)*(V(I + 9))/R(I) 
PRINT I; 

PRINT USING "##.###"; V(I); 
PRINT USING "MM# V(I + 9); 
PRINT USING "##.###"; P(I) 
NEXT I 


' Module average surface temperature. 


' Printing Voltages and Total Power to screen. 
' Calculating Power supplied by each Heater. 

' Voltage Drop across the Heaters. 

' Voltage Drop across the Precision Resistors. 


TOTPWR = 0! 

DO WHILE I < 10 'Loop for Total Power calculation. 

1 = 1+1 
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TOTPWR = TOTPWR + P(I - 1) 

LOOP 

PRINT "TOTAL POWER INPUT ="; 

PRINT USING "##.###"; TOTPWR 'Total PowerfWatts). 

NEXT I 


T3 = 0! 

PRINT TIMES; 'Printing the time to screen. 

DO WHILE I < 31 'Loop for calculating the average of 30 measurements 

taken for ' coolant inlet and outlet temperatures. 

1 = 1+1 

CHNLS = STR$(303) 

MSGS = "OUTPUT 9;MEAS TEMPT" + CHNLS 
PRINT #1, MSGS 
PRINT #1, "ENTER 9" 

INPUT #2, TMPS 
T(l) = VAL(TEMPS) 

T3 = T3 + T(I) 

TIN = T3 / 30 ' Average value for the coolant inlet temperature. 

LOOP 
NEXT I 
1 = 0 
T4 = 0! 

PRINT TIMES; 

DO WHILE I <31 
1 = 1+1 

CHNLS = STRS(304) 

MSGS = "OUTPUT 9;MEAS TEMPT" + CHNLS 
PRINT #1, MSGS 
PRINT #1, "ENTER 9" 

INPUT #2, TMPS 
T(l) = VAL(TEMPS) 

T4 = T4 + T(I) 

TOUT = T4 / 30 ' Average value for the coolant outlet temperature. 
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LOOP 

DIFF = TOUT - TIN ' Difference between inlet and outlet temperatures 

TAVG = (TIN + TOUT) / 2 'Average of inlet and outlet temperatures 
PRINT "Tin =", 

PRINT USING "MM", TIN; 

PRINT "Tout ="; 

PRINT USING "MM"-, TOUT; 

NEXT I 

FOR I = 0 TO 4 'Printing temperatures to screen. 

PRINT I; 

PRINT USING "M.M", TEMP(I), 

PRINT I + 5; 

PRINT USING "M.M"-, TEMP(I + 5); 

PRINT I + 10; 

PREMT USING "##,##"; TEMP(I + 10); 

PRINT I + 15; 

PRINT USING "M.M"-, TEMP(I + 15) 

NEXT I 

'Calculation for the heat transferred through the Insulation. 
DTINS = TEMP(l) - TEMP(O) 

K = .04 'Conductivity of the insulator 

A = .0219 'Surface area of the insulation Area. 

L = .0191 'Insulation Thickness. 

QINS = K * A * DTINS / L 

'Calculation for coolant thermophysical properties. 

A1 = -8.8216E-07 'Coefficients of density curvefit. 

B1 =-7.925 lE-04 
Cl =8.1055 

DENSITY = AI * TAVG ^ 2 + BI * TAVG + Cl 

A2 = 3.5284E-07 'Coefficients of specific heal curvefit.. 

B2 = .0010466 
C2 = .46955 

SPECHEAT = A2 * TAVG 2 + B2 * TAVG + C2 
QFLUID = DENSITY * SPECHEAT ♦ Flow * DIFF / 60 
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QTOTAL = QFLUID + QINS 
DELQ = QTOTAL - TOTPWR 
PRINT "Qfluid-'; 

PRINT USING ■'##.## "; QFLUID; 'Heat absorbed by the fluid. 

PRINT "Qins="; 

PRINT USING "##.## QINS; 

PRINT ''Qtotal="; 

PRINT USING "##.## "; QTOTAL; 

PRINT "DelQ="; 

PRINT USING "##.m "; DELQ 
PRINT "dTins="; 

PRINT USING "#.## "; DTINS; 

PRINT "dTfluid="; 

PRINT USING "##.m "; DTFLU; 

PRINT TIMES ' PRINT # I, "OUTPUT 9;BEEP" 'Make HP-3852 Beep when 

done 

PRINT # 1, "CLEAR 9" 'Clear and reset device 9. 

PRINT#!, "LOCAL9" 

INPUT " DO YOU WISH TO SAVE DATA ? ", SAVITS 'Decision to save the data 
or not. 

IF SAVITS = "n" THEN GOTO QUES 
IF SAVITS = "N" THEN GOTO QUES 


FILEIT: 'Storing the data into a file 

PRINT " IF A DATA FILE ALREADY EXISTS WITH THE NAME YOU SELECT, 

PRINT " PREVIOUS DATA WILL BE OVERWRITTEN " 

INPUT " ENTER THE NAME OF THE DATA FILE TO CREATE ", FILES 
IF FILES = "" THEN GOTO FILEIT OPENIT: 

OPEN FILES FOR OUTPUT AS #3 

PRINT #3, "mL/min" 'Flow rate , ml/min. 

PRINT #3, Flow 

PRINT #3, " Vhtr " 'Heater voltage drop. Volts. 

FOR I = 0 TO 8 

PRINT #3, USING "##,###"; V(I) 
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NEXT I 
PRINT #3, "Vr" 

Volts. 

FOR I = 0 TO 8 

PRINT #3, USING "##.###", V(I + 9) 
NEXT I 

PRINT #3, "POWER" 

FOR I = 0 TO 8 

PRINT #3, USING P(I) 

Watts. 

NEXT I 

PRINT #3, "TOT. POWER" 

PRINT #3, USING "##,###"; TOTPWR 
PRINT #3, "AVG SURF TEMP"; 
PRINT #3, USING "##.##"; TSURF 
temperature, ° C. 

PRINT #3, "Tin =" 

PRINT #3. USING "##.##"; TIN 
PRINT #3, "Tout =" 

PRINT #3, USING "##.##"; TOUT 
PRINT #3, "Qfluid=" 

PRINT #3, USING "##,##"; QFLUID 
coolant. Watts. 

PRINT m, "Qins=" 

PRINT #3, USING "##.##"; QINS 
Watts. 

PRINT #3, "Qtotal=" 

PRINT #3. USING "##.##"; QTOTAL 
PRINT #3, "DelQ=" 

PRINT #3, USING "##.##", DELQ 
PRINT #3, "dTins=" 

PRINT #3, USING DTINS 

insulation. 

PRINT #3, "dTfluid=" 


'Precision Resistor voltage drop, 

'Power supplied to each heater, 

'Total Power, Watts. 

'Module average surface 

'Coolant inlet temperature, "C. 
'Coolant outlet temperature, “C. 
'Heat absorbed by the 

'Heat loss through the insulation, 

'Total heat. Watts. 

'Difference in energy balance. Watts. 
'Temperature difference for the 
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'Temperature difference for the 


PRINT #3, USING DTFLU 

coolant. 

PRINT #3, "DEG C" 

FOR I = 0 TO 19 

PRINT #3, USING "##.##"; TEMP(I) 'Temperature readings, " C. 

PRINT #3, DATES 'Printing the date. 

PRINT #3, TIMES 'Printing the time. 

CLOSE #3 

QUES: 'Questioning. 

INPUT " DO YOU WISH TO TRY ANOTHER RUN ( Y/N ) ? ", AGAINS 
IF AGAINS = "Y" THEN GOTO STARTIT 
IF AGAINS = "y" THEN GOTO STARTIT 
IF AGAINS = "N" THEN END 
IF AGAINS = "n" THEN END 
GOTO QUES 
END 
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APPENDIX E. MATLAB CODES OF CALIBRATION CALCULATIONS 


The calibration calculations are performed by using Matlab codes as follows: 

% FLOWMETER (FISCHER & PORTER) CALIBRATION CURVE 
% 4ih ORDER POLYNOMIAL 

rate=[0 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 95 100]; % FLOW 
PERCENTAGE 

flow=[0 62.67 125.38 182.58 280.31 370,37 449.44 534.19 597.01 671.59 742.12 802.25 
884.17 941.18 995,02 1052.63 1120.45 1208.46 1267.43 \n4.Q(,]\%FLOW(ml/min) 
f=polyfit(rate,flow,4); 

% 4th ORDER POLYNOMIAL OF THE FLOWMETER CALIBRATION CURVE 

ratel=[0:l:100]; 

fl=polyval(f,ratel); 

plot(ratel,fl,rate,flow,'*'),grid; axis([-5 105 0 1500]); 
title ('FLOWMETER CALIBRATION CURVE'); 
xlabelCFLOW PERCENTAGE'); ylabel('FLOW RATE (mVmin)'); 
gtext('16 MAY 1995'); 

f; % COEFFICIENTS OF THE POLYNOMIAL 

% DENSITY CURVE OF THE COOLANT 
% BRAYCOMICRONIC889 
% MANUFACTURER'S DATA 

temp=[0 20 40 100 160]; % DEGREES CELSIUS 

dens=[.811 .794 111 .ITS .661]; % DENSITYg/ml 

f^polyfit(temp,dens,2); % 2nd ORDER POLYNOMIAL OF THE DENSITY 

CURVE temp 1=[0:1.0.160]; 

fl =polyval(f,temp 1), 

plot(temp 1 ,fl,temp,dens,"*'),grid; 

axis([-5 180 0.64 0.82]); 

title ('DENSITY CURVE of the COOLANT (BRAYCO MICRONIC 889)'); 
xlabelCTEMPERATURE (Deg. Celsius)'); ylabel(DENSITY (g/ml)'); 
f; % COEFFICIENTS OF THE POLYNOMIAL 
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% SPECIFIC HEA T vs TEMPERA PURE CURVE OF THE COOLANT 

% BRA YCO MICRONIC 889 

% DELSEN LAB'S DATA 

temp=[10 38 66 93]; % DEGREES CELSIUS 

heat=[.48 .51 .54 .57]; % SPECIFIC HEATcaUg-C 

f=polyfit(temp,heat,2); % 2nd ORDER POLYNOMIAL OF THE DENSITY 

CURVE temp l=[-20; 1.0; 100]; 

fl =polyval(f,temp 1), 

plot(temp 1 ,fl ,temp,heat,'*'),grid; 

axis([0 100 0.46 0.58]); 

title ('SPECIFIC HEAT CURVE'); 

xlabelCTEMPERATURE (Deg. Celsius)'); ylabel('SPECIFIC HEAT (cal/g-C)'); 
f; % COEFFICIENTS OF THE POLYNOMIAL 
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